clear;
clc;
H=19100e3;
Rz=6400e3;
Tobr=12*60*60;
Tmod=Tobr/2;
wobr=2*pi/Tobr;
c_light=3e8;
f_0=1602e6;
dt=5*60;
t=0:dt:Tmod;

rr=[0 Rz];

for i=1:length(t);
       rsv= (H+Rz)*[cos(wobr*t(i));sin(wobr*t(i))];
       plot([0 rsv(1)], [0 rsv(2)]);
       xlim([-H-Rz; H+Rz]);
       ylim([-Rz; H+Rz]);
       axis equal;
       pause(0.01);
       
       rv=rsv-rr;
end